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We discuss the solution to the Dark Matter problem provided by the 
lightest neutralino of the Minimal Supersymmetric Standard Model (MSSM) 
and highlight the role of the International Linear Collider (ILC) in deter- 



f- ^ , mining its cosmological relic density. 
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O ' 1. Introduction 

As deduced from the WMAP satellite measurement of the temperature 

anisotropies in the Cosmic Microwave Background, in combination with 

j^. data on the Hubble expansion and the density fluctuations in the universe, 

cold Dark Matter (DM) makes up ~ 25% of the energy of the universe pQ. 

The DM cosmological relic density is precisely measured to be 

^- ttvuh 2 = 0.113 ±0.009 (1) 



which leads to 0.087 < Odm h 2 < 0.138 at the 99% confidence level. In 
these equations, O = p/p c , where p c ~ 2 • 10~ 29 /i 2 g/cm 3 is the "critical" 
mass density that yields a flat universe, as favored by inflationary cosmology 
and as verified by the WMAP satellite itself; p < p c and p > p c correspond, 
respectively to an open and closed universe, i.e. a metric with negative 
or positive curvature. The dimensionless parameter h is the scaled Hubble 
constant describing the expansion of the universe. 

In the Minimal Supersymmetric Standard Model (MSSM), there is an 
ideal candidate for the weakly interacting massive particle (WIMP) which is 
expected to form this cold DM: the lightest neutralino Xi which is a mixture 
of the supersymmetric partners of the neutral gauge and Higgs bosons and 
is in general the Lightest Supersymmetric Particle (LSP). This electrically 
neutral particle is absolutely stable when a symmetry called R-parity is 
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conserved, is massive and thus non-relativistic or cold. Furthermore, it 
has only weak interactions and for a wide range of the MSSM parameter 
space, its annihilation rate into SM particles fulfills the requirement that the 
resulting cosmological relic density is within the range measured by WMAP. 
This is particularly the case in the widely studied minimal Supergravity 
(mSUGRA) scenario [2j and in some of its variants; see Ref. [3|. 

In this brief note, we discuss the contribution of the LSP neutralino 
to the overall matter density of the universe and highlight the role of the 
International Linear Collider (ILC) in determining this relic density. 

2. The Dark Matter relic density 

To derive the cosmological relic density, the standard treatment @] is 
based on the assumption [besides that the LSP should be effectively stable, 
i.e. its lifetime should be long compared to the age of the Universe, which 
holds in the MSSM with conserved R-parity that is discussed here] that 
the temperature of the Universe after the last period of entropy production 
must exceed ~ 10% of mo, an assumption which is quite natural in the 
framework of inflationary models 4 . In the early universe all particles were 
abundantly produced and were in thermal equilibrium through annihilation 
and production processes. The time evolution of the number density of the 
particles is governed by the Boltzmann equation 

dn x o 2 

— -i + 3Hn o =-(v cr ann )(n o - n^ ) (2) 

at A1 A i *i 

where v is the relative LSP velocity in their center-of-mass frame, cr a nn 
is the LSP annihilation cross section into SM particles and (. . .) denotes 
thermal averaging; no is the actual number density, while n e % is the thermal 

equilibrium number density. The Hubble term takes care of the decrease in 
number density due to the expansion, while the first and second terms on the 
right hand side represent, respectively, the decrease due to annihilation and 
the increase through creation by the inverse reactions. If the assumptions 
mentioned above hold, x? decouples from the thermal bath of SM particles 
at an inverse scaled temperature xp = m o/Tp which is given by [3] 

x F = 0.38M P (w ann )c(c + 2)m x o (g*x F y 1/2 (3) 

where Mp = 2.4- 10 18 GeV is the (reduced) Planck mass, g* the number of 
relativistic degrees of freedom which is typically g* ~ 80 at Tp, and c a 
numerical constant which is taken to be ^ ; one typically finds xp ~ 20 to 
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25. Today's LSP density in units of the critical density is then given by [3j 

^ = 2 -!i;f T /Ge T , w^ JM = r {va - ){x) dx ( 4) 



2 2.13 • 10 8 /GeV /• 00 (wa EL )(x) 

"v'i = — ; , 7 — 77 — r > Wlt h J(xp) = / tt 

x ^M P J(x F ) y ' J XF x 2 

Eqs. (JBlHQ provide an approximate solution of the Boltzmann equation 
which has been shown to describe the exact numerical solution very ac- 
curately for all known scenarios. Since Xi decouples at a temperature 
Tp <C m x , in most cases it is sufficient to use an expansion of the LSP 
annihilation rate in powers of the relative velocity between the LSPs 



va» 



= v cr{xiXi -► SM particles) = a + bv 2 + 0(v 4 ) (5) 



The entire dependence on the model parameters is then contained in the 
coefficients a and b, which essentially describe the LSP annihilation cross 
section from an initial S- and P-wave, since the expansion of the annihila- 
tion rate of eq. © is only up to 0(v 2 ). S-wave contributions start at 0(1) 
and contain 0(v 2 ) terms that contribute to eq. (jSJ) via interference with the 
0(1) terms. In contrast, P-wave matrix elements start at 0(v), so that 
only the leading term in the expansion is needed. There is no interference 
between S- and P-wave contributions, and hence no 0{v) terms. 

In generic scenarios the expansion eq. ® reproduces exact results quite 
well. However, it fails in some exceptional cases jS] all of which can be 
realized in some part of the MSSM parameter space, and even in mSUGRA: 

i) The expansion breaks down near the threshold for the production 
of heavy particles, where the cross section depends very sensitively on the 
cm. energy y/s. In particular, due to the non-vanishing kinetic energy of 
the neutralinos, annihilation into final states with mass exceeding twice the 
LSP mass ("sub-threshold annihilation") is possible. This is particularly 
important in the case of neutralino annihilation into W + W~ and hh pairs, 
for relatively light higgsino-like and mixed LSPs, respectively. 

ii) The expansion eq. (jSJ) also fails near s— channel poles, where the cross 
section again varies rapidly with y/s. In the MSSM, this happens if twice 
the LSP mass is near Mz, or near the mass of one of the neutral Higgs 
bosons. In models with universal gaugino masses, the Z- pole region is now 
excluded by chargino searches at LEP2 and we are left only with the Higgs 
pole regions which are important as will be seen later. 

Hi) If the mass splitting between the LSP and the next -to-lightest su- 
perparticle NLSP is less than a few times Tp, co-annihilation processes 
involving one LSP and one NLSP, or two NLSPs, can be important jSJ. 
Co-annihilation is important in three cases: higgsino or SU(2) gaugino like 
LSPs and when the LSP is degenerate in mass with fi or with the lightest 
top squark (the latter case hardly occurs in mSUGRA scenarios). 
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3. The relic density in the mSUGRA scenario 

The mSUGRA model [2] is the most widely studied implementation 
of the MSSM and it manages to describe phenomenologically acceptable 
spectra with only four parameters plus a sign: 

m , mi/2, A , tan/?, sign//. (6) 

where mo,mi/2 and Aq are the common soft SUSY-breaking terms of all 
scalar masses, gaugino masses and trilinear scalar interactions, defined at 
at the Grand Unification scale, tan (3 is the ratio of the vacuum expectation 
values (vev's) of the two Higgs doublets at the weak scale and \x is the 
supersymmetric higgs(ino) mass parameter. 

We use the Fortran code SUSPECT 6 to solve the RGE and to calculate 
the spectrum of physical sparticles and Higgs bosons, following the proce- 
dure outlined in Ref. T. In addition to leading to consistent electroweak 
symmetry breaking (EWSB), a given set of input parameters has to satisfy 
experimental constraints [SJ. The ones relevant for this study are [7]: 

- The total cross section for the production of any pair of sparticles at 
the highest LEP energy (209 GeV) must be less than 20 fb. 

- Searches for neutral Higgs bosons at LEP impose a lower bound on 
m^; allowing for a theoretical uncertainty, one requires m^ > 111 GeV. 

- Recent measurements of the muon magnetic moment lead to the con- 
straint on the SUSY contribution: —5.7 • 10 -10 < a^ SUSY < 4.7 • 10 -9 . 

- Allowing for experimental and theoretical errors, the branching ratio 
for radiative b decays should be 2.65 • 10~ 4 < B(b — ► sy) < 4.45 • 10 -4 . 

- Finally, the calculated Xi relic density has to be in the range Q. 

The output [7j is shown and partly commented in Fig.l. The black re- 
gions are those satisfying the DM constraint, eq. (1). In general, there are 
four familiar regions where this constraint is satisfied, i) Scenarios where 
both too and Vfiyi^ are rather small (the "bulk region") are most natural 
from the point of view of EWSB but are severely squeezed by lower bounds 
from searches for sparticles and Higgs bosons, ii) In the "co-annihilation" 
region one has to v o ~ rrif, , leading to enhanced destruction of sparticles 
since the f\ annihilation cross section is about ten times larger than that of 
the LSP; this requires m\/2 S> tuq. Hi) The "focus point" or "hyperbolical 
branch" region occurs at too 3> mj/2, and allows Xi to have a significant hig- 
gsino component, enhancing its annihilation cross sections into final states 
containing gauge and/or Higgs bosons; however, if nit is much higher than 
its current central value of 173 GeV, this solution requires multi-TeV scalar 
masses, iv) Finally, if tan/3 is large, the s— channel exchange of the CP-odd 
Higgs boson A can become nearly resonant, again leading to an acceptable 
relic density (the "A— pole" region). 
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Fig. 1. The mSUGRA {mx/2,mo) parameter space with all constraints imposed for 
Aq = 0,1-j. > and tan/3 = 10 (left) and 50 (right). The top quark mass is fixed to 
the new central value, m t — 172.7 GeV. The light grey regions are excluded by the 
requirement of correct electroweak symmetry breaking, or by sparticle search limits. 
In the dark grey regions f\ would be the LSP. The light pink regions are excluded by 
searches for neutral Higgs bosons at LEP, whereas the green regions are excluded 
by the b — ► 57 constraint. In the blue region, the SUSY contribution to g^ — 2 falls 
in the correct range whereas the red regions are compatible with having an SM-like 
Higgs boson near 115 GeV. Finally, the black regions satisfy the DM constraint. 

Recently, a fifth cosmologically acceptable region of mSUGRA param- 
eter space has been revived |0]. In a significant region of parameter space 
one has 2m < m^, so that s— channel h exchange is nearly resonant. This 
"h— pole" region featured prominently in early discussions of the DM den- 
sity in mSUGRA but seemed to be all but excluded by the combination of 
rising lower bounds on rrih and m from searches at LEP [S]. However, in 
recent years improved calculations ^Uj of the mass of the light CP-even h 
boson have resurrected this possibility for top mass values close to or larger 
than 178 GeV (thus, this region does not appear in Fig. 1). 



4. SUSY particle masses and the ILC 

Bounds on physical masses might be a more meaningful way to show 
the possibilities of mSUGRA than the ubiquitous plots of allowed regions 
in the space of basic input parameters, Fig. 1, in which one always fixes the 
values of some other free parameters (e.g., Aq, tan/3, m t ). One obtains the 
least biased view of the allowed ranges of masses by simply scanning over 
the entire parameter set that is consistent with a given set of constraints. 
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(s)particle 


mass bounds [GeV] 
Set I Set II Set III 


xf 


50 
105 
136 


53 
105 
137 


[53, 61] 

[105, 122] 

[280, -] 


h 


99 

114 
128 


99 

114 
128 


[630, -] 

[114, 122] 

[246, -] 


9 

d~R 


374 
444 
102 


383 

444 
110 


[383, 482] 
[774, -] 
[110, -] 



Table 1. Sparticle mass bounds in mSUGRA obtained by scanning over the entire 
allowed parameter space, defined by m t € [171 GeV, 185 GeV], (m^ + mj 2 )/2 < 2 
TeV, the lower bounds on sparticle and Higgs masses from collider experiments, the 
constraint on b — > sj, simple 'CCB' constraints and a conservative interpretation 
of the constraint from g^ — 2 (essentially the overlap of the 2a regions using r 
decay and e + e~ collider data). Set II adds the DM constraint to the above set of 
constraints. Set III is like Set II, except that the scanned region has been artifically 
limited to the h—pole region, where m^a < rrih/2. Only lower bounds are listed for 
Sets I and II, while for Set III the allowed range is given; a dash (-) means that 
the upper bound is directly set by the upper bound on the average stop mass. 



Table 1 lists lower bounds on the masses of some new (s)particles in 
mSUGRA, first without (set I) and then with (set II) the DM constraint. 
The lower bounds on many new (s)particles simply coincide with the bounds 
established by collider experiments. This is true for the lighter chargino, 
stau and scalar Higgs states, and essentially also holds true for the lighter 
stop. The bounds on the masses of the gluino and third neutralino arc 
essentially the same as that in a more general MSSM, as long as gaugino 
mass unification is maintained. Clearly the DM constraint still allows some 
new (s)particles to be quite light. One should emphasize, however, that 
usually the lower bounds in the table cannot be saturated simultaneously. 
Nevertheless, the possibility of light sparticles even in this simplest of all 
potentially realistic SUSY models that allow WIMP Dark Matter should be 
quite encouraging to experiments. 

Set III shows these bounds (including the DM constraint) when one 
confines oneself to the h— pole region discussed at the end of the previous 
subsection. In this case there are significant upper bounds on the masses 
of all gauginos. The reason is that one needs 2m^o ~ m^ < 120 GeV here, 
leading in mSUGRA with the assumed universality of gaugino masses at the 
GUT scale, to relatively light charginos and neutralinos as well as gluinos. 
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5. The determination of the relic density at the ILC 

Thus, in many scenarios SUSY particles can be produced abundantly at 
the next generation of high-energy colliders, in particular at the LHC and 
the ILC. However, to determine the predicted WIMP relic density (see the 
flowchart shown in the left-hand side of Fig. 2), one must experimentally 
constrain all processes contributing to the LSP pair annihilation cross sec- 
tion; this requires detailed knowledge not only of the LSP properties, but 
also of all other particles contributing to their annihilation. This is not a 
simple task and all unknown parameters entering the determination of Q x h 2 
need to be experimentally measured or shown to have marginal effects. 

Many high precision measurements are possible at the LHC ^J, but 
many of them can be vastly improved at the ILC |12l I13j . Because of the 
clean environment and the knowledge of the cm. energy of the initial beams, 
sparticle masses can be determined with high accuracy through kinematic 
endpoints and threshold scans. The results of one study ^3] i n a given 
mSUGRA scenario (SPSla point [H|) are summarized in the right-hand 
side of Fig. 2, where the achievable precision at collider experiments are 
compared with the satellite determination of Q x h 2 . The figure shows that 
the ILC will provide a part per mille determination of Q x h 2 in the case 
under study, matching WMAP and even the huge accuracy expected from 
Planck. The many possible implications of such measurements are outlined 
in the flowchart in the right-hand side of the figure. 
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Fig. 2. Left: Flowchart illustrating the possible implications of comparing Qhep> 
the predicted DM density determined from high energy physics, and Q C osmo, 
the actual DM density determined by WMAP and Planck. Right: Constraints 
in the (m x o, A(Q x h 2 )/Q x h 2 ) plane from the ILC and LHC. Constraints on 
A(Q x h 2 )/Q x h 2 from WMAP and Planck (which provide no constraints on m x o) 
are also shown; from Ref. JlSf . 
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Thus, if DM is composed of the lightest neutralinos, the LHC and par- 
ticularly the ILC will be able to determine the WIMP's properties and pin 
down its relic density. If these determinations match cosmological obser- 
vations to high precision, then (and only then) we will be able to claim to 
have determined what dark matter is. Such an achievement would be a 
great success of the particle physics/cosmology connection and would give 
us confidence in our understanding of the Universe. 
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